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ABSTRACT 

We present rest-frame optical images and spectra of the gravitationally lensed, star-forming galaxy 
J0900-I-2234 {z = 2.03). The observations were performed with the newly commissioned LUCIFERl 
near- infrared (NIR) instrument mounted on the Large Binocular Telescope (LBT). We fitted lens 
models to the rest-frame optical images and found the galaxy has an intrinsic effective radius of 
7.4 ± 0.8 kpc with a lens magnification factor of about 5 for the A and B components. We also 
discovered a new arc belonging to another lensed high-z source galaxy, which makes this lens system a 
potential double Einstein ring system. Using the high S/N rest-frame optical spectra covering H+K 
band, we detected H/3, [O III], Ha, [N II], and [S II] emission lines. Detailed physical properties 
of this high-z galaxy were derived. The extinction towards the ionized H II regions {Eg{B — V)) 
was computed from the flux ratio of Ha and H/3 and appears to be much higher than that towards 
stellar continuum {Es{B — V)), derived from the optical and NIR broad band photometry fitting. The 
metallicity was estimated using N2 and 03N2 indices. It is in the range of — ^ solar abundance, 
which is much lower than the typical z ^ 2 star-forming galaxies. From the flux ratio of [S II]A6717 
and [S II]A6732, we found that the electron number density of the H II regions in the high-z galaxy 
were ~ 1000 cm""^, consistent with other z ~ 2 galaxies but much higher than that in local H II 
regions. The star-formation rate was estimated via the Ha luminosity, after correction for the lens 
magnification, to be about 365 ± 69 Mq yr~^. Combining the FWHM of Ha emission lines and the 
half-light radius, we found the dynamical mass of the lensed galaxy is (5.8 ±0.9) x 10^" Mq. The gas 
mass is (5.1 ± 1.1) x 10^" Mq from the Ha flux surface density by using global Kennicutt-Schmidt 
Law, indicating a very high gas fraction of 0.79 ± 0.19 in J0900-I-2234. 
Subject headings: gravitational lensing — galaxies:high-redshift — galaxies: ISM 



1. INTRODUCTION 

The redshifts between 1 < z < 3 mark the time of peak 
cosmic star formation (Madau et al. 1998) and quasar 
activity (Fan et al. 2001). The current-day Hubble se- 
quence was being built up, and a large fraction of present- 
day stars were formed (Dickinson et al. 2003). Observa- 
tions of galaxies in this rcdshift range provide a direct 
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picture of the formation and evolution of galaxies and 
the assembly of their central supermassive black holes. 

By obtaining rest-frame optical spectra, we can per- 
form detailed studies of the star formation history, as 
well as properties of the ISM and stellar populations in 
galaxies. However, well-known diagnostic optical emis- 
sion lines are shifted into the near infrared (NIR) for 
z = 1 — 3. Since high-z galaxies (e.g., Lyman Break 
Galaxies or LBGs) are faint, it is difficult to obtain high 
signal to noise ratio (S/N) NIR spectra with current fa- 
cilities, especially for weak emission and absorption lines. 
Therefore, many previous NIR spectroscopic studies have 
been focused on the strong [O III] or Ha emission lines 
(e.g., Pettini et al. 2001; Erb et al. 2003, 2006b,c), from 
which only limited inferences can be drawn. 

One method of addressing this problem is to observe 
high-z star-forming galaxies which have been gravita- 
tionally lensed by foreground massive galaxies or clus- 
ters. Lensing can boost the observed flux of high redshift 
galaxies by several tens or more, so that high S/N spec- 
troscopy in the NIR becomes feasible. One of the first 
lensed galaxies to be studied in this way was the serendip- 
itously discovered MS1513-cB58 at z = 2.73 (Yee et al. 
1996; Bechtold et al. 1997; Teplitz et al. 2000). Sub- 
sequently, several groups began systematic searches for 
strongly lensed high redshift galaxies towards clusters of 
galaxies (Sand et al. 2005; Richard et al. 2008) and red 
galaxies in the Sloan Digital Sky Survey (SDSS, York et 
al. 2000) images (e.g., Belokurov et al. 2009; Kubo et 
al. 2009). These searches are finding lensed star- forming 



2 



Bian et al. 



galaxies at redshifts z > 2 (e.g., AUam et al. 2007; Lin 
et al. 2009; Diehl et al. 2009; Koester et al. 2010) wliieh 
can then be foUowed-up in the NIR, or rest-frame opti- 
cal and UV (e.g., Smail et al. 2007; Hainline et al. 2009; 
Finkelstein et al. 2009; Quider et al. 2009, 2010; Pettini 
et al. 2010). 

The object of this paper is the lensed system 
J0900-F2234 (z=2.03; Diehl et al. 2009) which was discov- 
ered in SDSS images, and then confirmed with follow-up 
deep optical g, r, i imaging and spectroscopy. The lens 
galaxy is at 2; = 0.49 based on SDSS spectroscopy. In 
this paper, we report J, H , Ks imaging and H+K spec- 
troscopy of J0900-|~2234 with the newly commissioned 
LUCIFER instrument (Mandel et al. 2008; Agcorges et 
al. in preparation 2010) at the Large Binocular Telescope 
(LBT; Hill et al. 2008). We detected the nebular emis- 
sion lines H^, [O III]AA4959,5007, Hq, [N II]A6583, and 
[S II]AA6717,6732 with simultaneous H+K band spec- 
tral coverage from 1.40-2.20/um. The wide simultaneous 
wavelength range coverage not only improves observing 
efficiency, but also helps to reduce the measurement un- 
certainties of the line ratios (e.g., Ha/Li/?) resulting from 
the flux calibration and slit losses variation between dif- 
ferent exposures. 

This paper is organized as follows: In § 2, we describe 
the observations and reduction of the LUCIFER data. 
In § 3, we describe the photometry, reconstruction of 
the source with lens modeling and line flux measure- 
ments. In § 4, we construct a simple lensing model us- 
ing the J0900-I-2234 NIR imaging. In § 5, we report 
the detailed physical properties of J0900-I-2234 and dis- 
cuss the systematic uncertainties of the physical prop- 
erty measurements. In § 6, our main results are summa- 
rized. Throughout this paper, we use a cosmology with 
Ho = 70 km s-i Mpc'S = 0.3, and = 0.7. The 
magnitudes are all AB magnitudes. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. LUCIFER 

LUCIFERl was built by a collaboration of five Ger- 
man institutes. It is the first of a pair of NIR im- 
agers/spectrographs for the LBT. It is mounted on the 
bent Gregorian focus of the left primary mirror. The 
wavelength coverage is from 0.85 to 2.4 fim {zJHK 
bands) in imaging, long-slit and multiobject spectroscopy 
modes . The detector is a RockweU HAWAn-2 Hd-CdTe 
2048 X 2048 pixeP array with a pixel size of 18.5 //m. 
The quantum efficiency (QE) is about 75% in the H and 
K bands. In the currently available seeing limited mode, 
the pixel scales arc 0.25" and 0.12" for N1.8 & N3.75 
cameras, respectively. The commissioning finished in 
early December 2009 and has been immediately followed 
by 12 nights of scientific demonstration observing. The 
instrument has been in regular science operation since 
mid-December 2009. 

2.2. LUCIFER Imaging & Spectroscopy 

J0900+2234 was observed with LUCIFERl on the 
LBT on 6 January 2010 in imaging mode and on 13 
February 2010 in long slit spectroscopy mode (Table 1). 
In the imaging mode, the N3.75 camera was used; the 
18.5 lira pixels correspond to 0.12" on the sky, for a 
4.0' X 4.0' field-of-view. The seeing was 0.7" with thin 



clouds. For each filter (J, H, and Ks), we obtained ten 
randomly dithered exposures of 60s each. For the spec- 
troscopic observations, the seeing was 0.5-0.6" and the 
conditions were photometric. The spectra were taken 
with N1.8 camera, yielding a plate scale of 0.25" per 
pixel. An order separation filter and the 200 1/mm H+K 
grating were used to cover the wavelength range from 
1.40 to 2.20 iJLva. simultaneously. A 1.0" by 3.9' slit was 
used, resulting in spectral resolution of 16 A. The total 
exposure time was 24 x 300 seconds with the telescope 
nodded ~ 7" along the slit between each exposure. The 
instrument was rotated to PA=86.57° in order to place 
the two brightest lensed components in the slit. Dark 
exposures and internal quartz lamp fiats were obtained 
in the afternoon, and an Ar lamp was observed for wave- 
length calibration during the night. 

2.3. Data Reduction 
2.3.1. Imaging Reduction 

The imaging data were reduced using standard IDL 
routines. The dark frames were median combined to 
create a master dark frame. A super-sky fiat was con- 
structed from combining the science frames after dark 
subtraction. The science images were dark subtracted, 
and then divided by the super-sky fiat to correct for the 
detector response. The science images were registered 
to the SDSS-DR6 catalog using SCAMP (Bertin 2006) 
for astrometric calibration; the residuals are less than 
0.1" RMS. Finally, the sky background was subtracted 
from the science frames which were then co-added using 
SWarp (Bertin et al. 2002). Photometric zeropoints were 
determined from 2MASS stars in the field-of-view. 

The fully reduced J, H, Ks-hand images arc shown in 
Fig. 1. In addition to the lensed galaxy (components A, 
B, C, and D in Fig. 1) found by Diehl et al. (2009), we 
discovered another component, Al. The color of Al sug- 
gests that it is also a high-redshift galaxy image lensed 
by the foreground red galaxy, but not the same galaxy 
appearing in images as A, B, C or D (see details in sec- 
tion 3.1). 

2.3.2. Spectroscopy Reduction 

To reduce the LUCIFER spectra, we modified an IDL 
long-slit reduction package written by G. D. Becker for 
NIRSPEC (Becker et al. 2009). In this package, the fol- 
lowing sky background subtraction procedure is used. 
First, four types of calibration files were created: a 
median-combined normalized internal flat field to correct 
the pixel-to-pixel variations in QE, a median-combined 
dark image, and two transformation maps created by 
tracing bright standard stars and bright sky lines. These 
two transformation maps were used to transform the NIR 
detector (x, y) coordinates to the slit position and wave- 
length coordinates. Then, the dark image was subtracted 
from the science frames which were then divided by the 
fiat field. The background in each science frame was sub- 
tracted using a fit to the science frames obtained before 
and after it in two dimensions, based on the transforma- 
tion maps. A second-order polynomial function and a b- 
spline function were used to fit the sky background resid- 
ual along the slit direction and the dispersion direction, 
respectively. After the sky background subtraction, the x 
and y direction distortion of the science frames were cor- 
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rected based on the transformation maps created in the 
first step. The one-dimensional spectra were extracted 
from individual two-dimensional spectrum frames and 
combined together to give the averaged spectrum shown 
in Fig. 2. The wavelength calibration was derived from 
observations of Ar lamps and applied to the averaged 
spectrum. The spectrum was corrected for telluric fea- 
tures with the spectrum of an AV star observed at the 
same average air mass as the science exposures. Finally, 
we derived the flux calibration by normalizing the spec- 
trum to the if-band magnitude. 

3. MEASUREMENTS 

3.1. NIR Photometry: Discovery of a new Lensed 

Galaxy 

The GALFIT program (Peng et al. 2002, 2010) was 
used to model the central lensing galaxies and lensed 
components. We fit each of the two central massive 
galaxies with de Vaucouleurs models, and fit the three 
lensed knots and two arcs with exponential disk models. 
The initial input parameters - namely, position, total 
magnitude, axis ratio, and position angle for each com- 
ponent - were determined using SExtractor (Bertin & 
Arnouts 1996). We used the Ks band co-added image 
to fit the profile parameters (Fig. 3 and Table 2), and 
fixed these parameters when fitting the J and _ff-band 
images to obtain the magnitude of each component (Ta- 
ble 3). The values of best-fit reduced for J, H, and 
Ks bands are 1.064, 1.094, and 1.088, respectively. 

The NIR colors of the A and B components are similar, 
with J — _ff ~ 0.1 and H — Ks ^ 0.0, while the colors 
of the C and D components appear to be different from 
those of the A and B components. For D component, the 
reason for this discrepancy may be that the exponential 
disk profile is not a good fit to the observed arc struc- 
ture. For the C component, the follow-up spectroscopy 
was performed with MMT/Blue Channel spectrograph, 
which shows that the C component is not part of this 
lens system. 

Current data are not yet deep enough to uncover other 
components of the weak Al component. If confirmed by 

spectroscopy and deeper imaging, this newly discovered 
arc would make J0900-I-2234 a rare example of a double 
Einstein ring system lensed by galaxy. 

3.2. Emission Line Measurements 

We used the IDL MFFIT^" package to fit the the emis- 
sion lines. The results are listed in Table 4. The H/3, 
[O III]A4959 and [O III]A5007 were fit as Gaussian func- 
tions individually. The positions of Ha and [N II]A6583 
as well as [S II]AA6717, 6732 are close to each other, so 
they were fit by two Gaussian functions together to de- 
blend these two lines. The observed line wavelength, line 
flux, and full width at half maximum (FWHM) were de- 
rived from the Gaussian fitting (Table 4). Monte Carlo 
(MC) simulation was used to estimate the uncertainties 
of line flux and FWHM. One thousand artiflcial spectra 
were generated by perturbing the flux of each data point 
from the true spectrum by a random amount propor- 
tional to the la flux error. We then measured the line 
flux of each fake spectrum with leaving all the parame- 
ters free for both the single lines and deblened line pairs. 

http://cow.physics.wisc.edu/ craigm/idl/idl.html 



The standard deviations of the distributions of the line 
flux and the FWHM were adopted as the uncertainty 
of measurements of the line flux and the FWHM (Ta- 
ble 4). The emission line redshifts of A and B knots are 
2.0321±0.0009 and 2.0318±0.0005, respectively, and are 
consistent with each other within Az/{1 + z) = 0.0001. 
These results also agree with the mean redshift from the 
optical spectra, z = 2.0325 ± 0.0003 (Diehl et al. 2009). 

4. LENSED MODEL 

The lens modeling is performed using LENSFIT (Peng 
et al. 2006). Briefly, LENSFIT is patterned after, and 
works like GALFIT (Peng et al. 2002, 2010). LENSFIT 
is designed to specifically deal with situations where the 
image geometry is crowded: it allows one to decompose 
foreground and background galaxy light profiles, and de- 
termine the lens deflection model, simultaneously. The 
number of light proflle and singular isothermal ellipsoid 
(SIE) deflection models is unrestricted, and the optimiza- 
tion process is done using the Levenberg-Marquardt al- 
gorithm in Numerical Recipes (Press et al. 1992). 

In J0900-I-2234, the deflection model is potentially 
quite complicated because the system is embedded in a 
compact cluster environment. There are at least two pri- 
mary deflectors, and potentially 5 in all, just within the 
Einstein ring of the system which has a radius of roughly 
7.7 arcseconds. Due to the low signal-to-noisc of the 
source detection, it is difficult to determine an accurate 
lens model. However, it is clear from visual inspection 
that a single SIE deflector probably would not suffice be- 
cause the geometric center of the lensed arcs falls in the 
gap between the two primary deflectors. Therefore, the 
simplest lens deflection model we adopt is that of two 
SIEs held to the position of the light profile of the two 
primary lensing galaxies. We also hold the SIE axis ratios 
fixed to the light profile models of the lenses, but allow 
the Einstein ring radii and position angle parameters to 
be free in the fit. We simultaneously optimize the light 
profiles of all the foreground sources and the background 
galaxy light profiles, in all 14 objects in J-band image. 

With the limitation of our lens model and the data S /N 
in mind, we infer that the size of the background source is 
0.8-1.0 arcsecond which corresponds to the physical size 
of 6.6-8.2 kpc, in deprojected angular size, and it has a 
luminosity of 21.7 to 21.8 mag (AB) in J band, given a 
magnification of the A and B components together is 4.6- 
5.0. For the further analysis, we use the effective radius 
of 7.4 kpc and the A and B component magnification of 
4.8. The Sersic index is quite high, from n=3.5 to n=4. 
The major uncertainty of the lens fitting is caused by 
the influence of the bright foreground source C, which is 
right on the Einstein ring and the source elsewhere has 
quite low S/N. 

5. PHYSICAL PROPERTIES 

With LUCIFER we obtained H+K (1.40-2.20 /xm) 
spectra of the A and B components of the lensed galaxy 
J0900-F2234 (Fig. 2). The H/3, [O III]AA4959,5007, Ha, 
[N II]A6583, and [S II]AA6717,6732 lines were detected. 
In this section, we discuss the physical properties which 
can be derived from these observations. 

5.1. Emission-Line diagnostics 
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First, it is important to examine the observed values 
of well-known empirical line ratio diagnostics to deter- 
mine whether the emission is from an H II region, active 
galaxy, or shocked gas. The line ratios of the A and B 
components on the diagram of Baldwin, Phillips, & Ter- 
levich (1981) (BPT diagram) are shown in Fig. 5. The 
weak N II emission rules out the possibility that the emis- 
sion is from the activity of a central active galactic nu- 
cleus (AGN). The location of the A and B components on 
the BPT diagram are similar to other high-z star form- 
ing galaxies (e.g., Liu et al. 2008; Hainline et al. 2009). 
They have relatively higher [O III]A5007/H/3 values and 
lower [N IIA6584/Hq; values compared to the local star- 
forming galaxies from SDSS (Fig. 5). This offset could be 
the result of relatively intensive star formation activity 
and high electron density in these high redshift galax- 
ies, compared to the local sample. Overall, we conclude 
that the line emission is from gas that is photo-ionized 
by a hot stellar continuum, that is, the emission is from 
star-formation regions. 

5.2. Reddening and Extinction 

We estimated the extinction of the lensed galaxy by 
two methods: (1) a fit to the stellar continuum as mea- 
sured by broad band photometry, and (2) Balmer decre- 
ment. These two methods represent dust extinction to 
the stellar continuum Eg {B — V) and the dust extinction 
to the nebular gas in the H II region {Eg{B — V)), re- 
spectively. One expects these two estimates to differ, and 
indeed some previous studies have observed that there is 
more extinction towards the ionized gas than the stellar 
continuum in local star-forming galaxies (Calzetti et al. 
2000) and high redshift galaxies (e.g., Forster Schreiber 
et al. 2009). In contrast, Erb et al. (2006b) did not find 
any extinction difference between ionized gas and the 
stellar continuum in star-forming galaxies at redshift of 
~ 2. For lensed galaxies, the reddening derived from 
individual images often differs because of reddening by 
different amounts of foreground dust, since the images 
traverse different sight-lines through the lensing galax- 
ies. 

With optical photometry alone, one can measure only 
the UV continuum shape in high redshift galaxies, mak- 
ing it difficult to distinguish between reddening by dust 
and age of the stellar population. By adding the NIR 
photometry, we can measure the Balmer break (3600- 
3700 A) which is sensitive to the age of the stellar pop- 
ulation, and the 4000 A break strength only weakly de- 
pends on metallicity. By combining our NIR photometry 
with the optical, we can fit stellar population models and 
simultaneously determine the age of the galaxy and av- 
erage extinction. 

Here we used the J, H, and Ks-hand (rest-frame op- 
tical bands) magnitudes combining with the g, r, and i- 
band (rest-frame UV bands) magnitudes (3" from Diehl 
et al. (2009) to estimate the Es{B - V) and the age of 
the galaxy The Bruzual & Chariot (2003, BC03) stan- 
dard simple stellar population (SSP) model was used to 
build the spectral templates of star-forming galaxies. We 
adopted constant star formation rate (SFR) with the 
Chabrier initial mass function (IMF) Chabrier (2003), 
and solar metallicity, to generate a series of spectra with 
different ages and reddening. Fig. 4 shows the best-fit 



spectra with the photometry data of A and B compo- 
nents. From the fitting, we found that the ages of the 
components A and B were consistent with each other, 
and equal to 180 Myr, with the values of Es{B — V) 
equal to 0.07 and 0.20 for the components A and B, re- 
spectively. Based on the metallicity found in section 5.4, 
we also generated the spectra with 0.25 solar metallicity 
to fit the broad band photometry data, and found similar 
results. 

With this fit to the stellar population we derive the 
intrinsic stellar mass of the source galaxy is 1.9 x 
lO^'' Mq after correcting for lensing magnification, which 
is smaller than the mean stellar mass of 3.6 x lO^'^ Mq 
found in a z 2 UV-selected galaxy sample (Erb et 
al. 2006c). Using Ks band magnitude to approach the 
rest-frame R band magnitude, we found the M/Lr in 
J0900-h2234 is 0.13. Shapley et al. (2005) found that the 
mean of M/Ln is 0.29 with a large scatter from 0.02 to 
1.4 (where the values have been divided by 1.8 to con- 
vert the Salpeter IMF to Chabrier IMF) in z - 2 UV- 
selected star-forming galaxies. Considering the high scat- 
ter, M/Lr in J0900-h2234 agrees with other z - 2 star- 
forming galaxies. The intrinsic UV (~ 1700 A) bright- 
ness of the J0900-I-2234 is 22.1 magnitude, so that the ab- 
solute magnitude is M(1700)ab = -22.7, compared to L* 
= -20.20—20.97 at z ~ 2 (e.g. Oesch et al. 2010; Reddy 
et al. 2008; Reddy & Steidel 2009). Thus, J0900-K2234 
has L/L* « 5—10, and is an intrinsically luminous galaxy. 
Galaxies with these high luminosity is very rare, and 
their number density is about a few 10~^ Mpc~^ mag~^ 
base on the luminosity function from Reddy et al. (2008) 
and Reddy & Steidel (2009). The high intrinsic UV lu- 
minsity indicates a high SFR (see details in Section 5.3) 
in this galaxy. 

We estimated Eg{B - V) of J0900-I-2234 using the fiux 
ratio of Ha and H/3 lines. Under Case B H I recombina- 
tion at a temperature T=10,000 K and electron density 
n <100 cm"'^ (Zaritsky et al. 1994), the intrinsic value 
of Ha/H/3 is 2.86 (Osterbrock & Ferland 2006). Using 
the extinction law proposed by Calzetti et al. (2000), 
we found the values of Eg{B - V) were 0.84 ± 0.31 and 
0.59 ± 0.08 for the A and B components. For both A 
and B components, the values of Eg{B — V) are signifi- 
cantly larger than those of Es{B — V), in agreement with 
the results of Calzetti et al. (2000) and Forster Schreiber 
et al. (2009). The stars and line-emitting gas are not 
co-spatial, and the line-emitting regions are dustier than 
the galaxy as a whole. Note that the low S/N of the H/3 
line in knot a makes the uncertainty of the Eg{B — V) 
for component A large. Therefore, we use the value of 
Eg{B -V)= 0.59 ± 0.08 of knot B for further analysis. 

5.3. Star Formation Rate 

We estimated the star-formation rate (SFR) in 
J0900-I-2234 using the luminosity of Ha (I/Ha, Kennicutt 
et al. 1998). The relation between the SFR and Lho is 

SFR(MQyr-i) = 7.9 x io-42 Lh« 

ergs s ^ 

The Ha fluxes of the A and B components are 
37.01 ± 0.79 X 10^1^ ergs s'^ cm'^ and 84.64 ± 3.21 x 
10-1^ ergs s-i cm-2. We used EgiB^-V) = 0.59±0.08to 
correct the fluxes for extinction and found the extinction- 
corrected fluxes are (225.5 ±55.5) x 10~^^ ergs s~^ cm~^ 
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and (516.0 ± 128.0) x lO'^"^ ergs s^^ cni"2_ ^he 
extinction-corrected luminosities are therefore (15.4 ± 
3.8) X 10''3 ergs s'^ and (6.75 ± 1.66) x lO''^ ergs s'^. 
The lensed SFRs of components A and B are 533 ± 
131 Mq yr-i and 1220 ± 303 Mq yr"!. We added the 
components A and B SFRs together, corrected for the 
magnification, and found that the intrinsic SFR of the 
lensed galaxy was 365 ± 69 Mq yr~^. Note that if we 
use the Chabrier IMF, the SFR derived from the Lh^ 
will be 203 ± 38 Mq yr"!. We also estimate the SFR 
from the UV luminosity (Kennicutt et al. 1998) which 
is derived from the average of g and r band flux. We 
found that the SFR is 276 Mq yr for Shalpeter IMF and 
153 Mq yr for Chabrier IMF by correcting the dust ex- 
tinction, Es{B -V) = 0.07(0.20) for A (B) components. 
The SFR based on the dust-corrected Ha emission is con- 
sistent with that from the dust-corrected UV emission. 
If we use Es{B — V)= 0.20 to correct the extinction of 
Ha flux, the Ha-derived SFR of 110 Mq yr"! will be 
significant lower than UV-derived SFR. This supports 
the result that the Eg{B - V) is larger than Es{B - V). 
The average of the SFR for z ~ 2 star-forming galaxies is 
(SFRHa)= 31 Mq yr-i and (SFRuv)= 29 Mq yr-i(Erb 
et al. 2006c). Despite the uncertainties in interpreting 
the SFR measurements, we conclude that the SFR of 
J0900+2234 is an order of magnitude higher than typi- 
cal star-forming galaxies at 2; ~ 2. 

5.4. Oxygen Abundance 

Emission lines from H II regions can be used to mea- 
sure metallicity. The oxygen abundance was calcu- 
lated using the following indicators: the N2 index (N2= 
log(i^[Nii]A6584/i^Ha)) and the 03N2 index (03N2= 

log{(-F[oiii]A5007/-FH/3)/(i^[Nii]A6584/^Ha)}), which have 
been well calibrated using nearby extragalactic H II re- 
gions to measure 0/H (Pettini & Pagel 2004). The ad- 
vantage of using the N2 and 03N2 indices as metal alnm- 
dance estimators is that Ha and the [N II]A6584 pair, as 
well as the H^ and [O III]A5007 pair, are close in wave- 
length, making these two estimators relatively insensitive 
to extinction. 

The N2 index is sensitive to the oxygen abundance 
Storchi-Bergmann et al. (1994), and was further cali- 
brated by Raimann et al. (2000) and Denicolo et al. 
(2002). Here we use the best linear fit between the N2 
and (12 -hlog(0/H)) from Pettini & Pagel (2004), 

12 -h log(0/H) = 8.90 -h 0.57 X N2, (2) 

where N2 is in the range from -2.5 to -0.3, and the la 
error of the measurements of log(0/H) is 0.18. We find 
the values of 12 -|- log(0/H) equal to 8.12 ± 0.21 and 
8.23 ± 0.19 for components A and B, respectively, which 
are about 0.27 ± 0.13 to 0.35 ± 0.15 solar abundance (for 
the sun: 12 + log(O/H0) = 8.69 (Asplund et al. 2009)). 

The 03N2 index was first introduced by Alloin et al. 
(1979), and was well calibrated by Pettini & Pagel (2004) 
using 137 extragalactic H II regions. The relation be- 
tween the metallicity (12 -Mog(0/H)) and the 03N2 in- 
dex is: 

12 -I- log(0/H) = 8.73 - 0.32 x 03N2, (3) 

where 03N2 is in the range from -1 to 1.9, and the Icr 
error of the measurements of log(0/H) is 0.14, which 
indicates a less scattering than the N2 index. We found 



the values of 12 + log (0/H) of components A and B were 
8.00 ± 0.16 and 8.09 ± 0.15, which are 0.21 ± 0.08 and 
0.25 ± 0.08 solar abundance. 

Both the N2 and 03N2 indices indicate that 
J0900+2234 has a low oxygen abundance compared with 
other 2 ~ 2 star-forming galaxies. Fig. 6 shows the stel- 
lar mass and metallicity relation in the local SDSS star- 
burst galaxies and z ^ 2 star-forming galaxies (Erb et 
al. 2006a). The metallicity is derived from N2 index. 
We averaged the value 12 -|- log(0/H) of the A and B 
components and the stellar mass is from the broad band 
fitting result. Our data point is significantly lower than 
the mass and metallicity relation in star-forming galaxies 
at 2; ~ 2(Erb et al. 2006a). Mannucci et al. (2010) pro- 
posed a more general fundamental relation between stel- 
lar mass, SFR and metallicity derived with SDSS galax- 
ies, in which the metallicity decreases with the increase 
of SFR. We used the relation derived in Mannucci et 
al. (2010) (equation 2) to obtain the predicted metallic- 
ity from SFR and stellar mass and found the value is 
8.45. Although our measured metallicity is ~ 0.27 dex 
lower than the predicted value, this dispersion is con- 
sistent with the results of Mannucci et al. (2010) that 
the distant galaxies show 0.2-0.3 dex dispersions for the 
fundamental metallicity relation. 

5.5. Electron Density 

The flux ratio of [S II]A6717 and [S II]A6732 
was used to estimate the electron density in the 
H II region of J0900+2234. The IRAF task 
stdas. analysis. nebular. temden was used to compute 
the density. We assumed that the temperature in 
these regions were 10.000 K. We found the values of 
^[sii]A67i7/J^[sii]A6732 Were 0.88 ± 0.24 and 0.84 ± 0.31 
for the components A and B. The large error bar is due 
to the low S/N (< 5) of these [S II] lines in both aper- 
tures. The values of the ratio yielded a range of electron 
number density lQ2dtf^ cm'^ and 1166l^f cm'^, re- 
spectively for these two components. Nonetheless, these 
densities are much higher than those typical of local H II 
regions ^IQO cmT^ (e.g., Zaritsky et al. 1994). The 
electron density of J0900-I-2234 is similar to that in the 
lensed galaxies, the Cosmic Horseshoe and the Clone, 
which is also derived from the ratio of [S II] double 
lines (Hainline et al. 2009). The high electron density 
was also found in the Cosmic Horseshoe with 5 000—25 
000 cmr^ derived from the ratio of C III]AA1906, 1908 
doublet (Quider et al. 2009). The high electron density 
implies the compact size of the H II regions in the high-z 
galaxies. If these high-z H II regions follow the similar 
electron density-size relation found in the local galaxies 
(Kim & Koo 2001), their sizes should be less than 1 pc. 

5.6. Virial Mass and Gas Mass 

The width of the emission lines can be used to probe 

the dynamics and the total mass of the parent galaxy. By 
fitting the emission lines (e.g.. Ha, [O III]) as Gaussian 
profiles, we found that the FWHMs of Ha of A and B 
knots were 19.0Aand21.8A. These observed FWHMs 
were corrected by subtracting the instrument resolu- 
tion 16 A) in quadrature. The corrected FWHMs 
were converted to the 1-D velocity dispersion with a = 
FWHM/2.355 x c/A. The values of a for the A and 
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B knots arc 66 ± 5 km s^^ and 95 ± 6 km s^^ respec- 
tively. We averaged the a components of A and B to 
estimate the velocity dispersion of J0900-t-2234, which is 
81 ±4 km s~^. Our velocity dispersion is about 30% lower 
than that in a sample of 2; ~ 2 star- forming galaxies (Erb 
et al. 2003), who found a mean velocity dispersion of 
- 110 km s~i. 

The mass of the galaxies can be estimated by assuming 
a simplified case of a uniform sphere (Pettini et al. 2001) 



M.. = 1.2xlO-Me(^^5^^)^^ (5) 

where G is the gravitational constant and ri/2 is the 
half-light radius, which is 7.4 ± 0.8 kpc from the lens 
model. We find that the dynamical mass of J0900-t-2234 
is (5.8 ± 0.9) X 10^° Mq, an upper limit since the line 
emitting gas may not be in Virial equilibrium with the 
gravitational potential of the galaxy. 

Using the global Kennicutt-Schmidt law (Kennicutt et 
al. 1998), we convert the surface gas density (Sgas) with 
the SFR derived from Ha: 

^.«« = l-6xlO--(^;^^)°'"Me.c- (6) 

where TiHa — LnalTt is the surface density of Ha lumi- 
nosity (Erb et al. 2006b; Finkelstein et al. 2009). Then 
the gas mass, Mgas- can be derived from 'Shu x r"^, 
which is 5.1 ± 1.1 x 10^" Mq. And the gas fraction, 

fgas = Mgas/{Mgas + M.tellar), is 0.74 ± 0.19. Erb Ct 

al. (2006a) found the gas fraction increases when the 
stellar mass decreases in the UV-selected star forming 
galaxies at z ~ 2. The gas fraction in J0900-I-2234 is 
higher than the mean gas fraction (0.48 ± 0.19) of Erb 
et al. (2006a) galaxies in the similar stellar mass bin 
(1.5 ± 0.3 X 10^° Mq), but still Icr consistent. The total 
baryonic mass Mgas + Mgteiar is 6.9 ± 1.1 x 10^° Mq, 
which is in the agreement with the dynamic mass. 

5.7. Measurement Uncertainty 

The main uncertainty of the measurement of the ex- 
tinction and abundance is from the uncertainties in the 
measurements of the weak emission lines (e.g. H/3 and 
S II) in the observed NIR. There are many prominent 
sky emission lines in the NIR, especially in the H band, 
which can contaminate our line measurements, especially 
for the weak emission lines. In the LUCIFER spec- 
tra of J0900-I-2234, the H/3 sits near the blue cut of 
the H band, where there are several moderate strong 
sky emission lines and the instrument efficiency, which 
contributes large flux errors (S/N< 5) and underesti- 
mate/overestimate the fluxes. The higher Eg{B — V) 
value of the A (fainter) component and systematically 
lower metallicities estimated from 03N2 index compared 
to those from N2 index may imply the flux of H/3 is un- 
derestimated. 

The systematic uncertainty of SFR, associated with 
the absolute flux of the Ha, is mainly due to the uncer- 
tainties in the absolute calibration of the NIR spectra. 
We calibrated the spectra with an AV star, and scaled the 
spectral flux to the H band magnitude to correct the light 



loss in the spectrograph slit. The scale factors are ^1.6 
(~ 1.8) for the A (B) knots, which is a reasonable value 
for light loss. Converting the lensed Ha luminosity to the 
intrinsic Ha luminosity is another uncertainty for SFR 
estimation. The uncertainty of the lens model fltting is 
20 — 30%. The uncertainty also affect the estimation of 
the stellar mass and dynamical mass, which depend on 
the results of the lens model. The Ha and SFR relation 
(Kennicutt et al. 1998) is derived from stellar synthesis 
models with solar metallicity. The SFR in J0900-I-2234 
could be overestimated, given the low metallicity (^ 0.25 
Zq) in this galaxy. 

6. CONCLUSIONS 

We present LBT/LUCIFERl NIR (rest-frame opti- 
cal) imaging and spectroscopy of the lensed galaxy 
J0900-h2234 (z = 2.03). The lensed components A and 
B were placed in the slit to obtain NIR spectra covered 
from 1.40 to 2.20 /xm. The detailed physical properties 
of the lensed star-forming galaxy were studied using the 
rest-frame optical spectra (Table 5). We summarize the 
main results as follows. 

• The new imaging was used to construct a lensing 

model. The magniflcation factor for the fluxes of 
the lensed galaxies (A plus B) is estimated to be 
4.8. 

• A new lensed arc Al was discovered by the deep 
NIR J, H, and i^s-band images. The colors and 
position imply that this arc does not have the same 
source galaxy as the other four components discov- 
ered before, suggesting the presence of a rare dou- 
ble Einstein ring. Follow-up spectroscopy shows 
that the C component does not have the same 
source galaxy as A and B. 

• We flt the optical and NIR broad band photometry 
to theoretical stellar population spectral templates 
(BC03), and found the galaxy age to be 180 Myr, 
EsiB -V) = 0.07 (EsiB -V) = 0.20) for the A 
(B) component. The stellar mass is 1.9x10^° Mq, 
and the intrinsic luminosity of the galaxy is L/L* 
w 5 - 10. 

• Using the flux ratio of Ha and H^, we found 

that the emission line gas suffers an extinction of 
Eg{B-V) is 0.84 ±0.31 (0.59 ±0.08) for the A (B) 
component, which is much higher than Es{B — V). 
This result implies that there is more extinction to- 
wards the ionized gas than the stellar continuum. 

• The oxygen abundance was estimated using the N2 

and 03N2 indices. We found that the metallic- 
ity in this star forming galaxy is only 21 — 35% 
solar, which is somewhat lower than typical star- 
forming galaxies at 2; ~ 2. The electron number 
density, Ue, was measured using the flux ratio of 
[S II]A6717 and [S II]A6732. The Ue is in the range 
of 10291^333 cni-3 (1166l™f pm-S) for the A (B) 
component. These values are consistent with other 
studies of high-z galaxies (e.g., Hainline et al. 2009; 
Quider et al. 2009), but much larger than that in 
local H II regions. 
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• The extinction-corrected luminosity of Ha emission 
line was used to estimate the SFR, and the intrin- 
sic SFR of J0900-I-2234 is 365 ± 69 Mq yr"! which 
is one magnitude higher than a typical z ~ 2.0 
UV-selected star-forming galaxy. We estimate the 
Virial mass of J0900-h2234 from the FWHM of the 
Ha emission line and find a dynamical mass of 
(5.8 ± 0.9) X 10^° Mq. The gas mass of the galaxy 
is (5.1 ±1.1) X 10^° Mq estimated from the Ha flux 
by assuming the global Kennicutt-Schmidt law. 

Comparison to the local star-forming galaxies, the 
J0900+2234 (z=2.03) features as much higher electron 
density. The high electron density and high ionization 
state (e.g., Hainline et al. 2009) probably cause the off- 
set between the H II regions in the high-z galaxies and 
those in the local galaxies in the BPT diagram. This 
star-forming galaxy also has a much higher SFR, and a 
relatively lower age, stellar mass, and metallicity than 
those in the typical high-z star-forming galaxies. 



The commissioning of LUCIFERl on the 8.4m LBT 
left primary mirror (SX) has allowed a detailed study of 
the z = 2 galaxy J0900-I-2234 in relatively short expo- 
sure times. We look forward to future NIR observations 
of other lensed galaxies, as well as field galaxies which 
are not lensed, which will give a similarly richly detailed 
picture of galaxies in this important stage of galaxy evo- 
lution. 
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Fig. 1.— J, H, and Ks band images of ,10900+2234 from LUCIFER, North is up, and cast is left. Each image shown has 40" X 40" 
The previously discovered A, B, C, and D components are labeled in the J-band image. The A, B, and D components are images of the 
star-froming galaxy at z ^ 2. Our follow-up spectroscopy shows that C is not an image of the same galaxy. The newly discovered Al arc 
from another lensed galaxy is also labeled in the image. 
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Fig. 2. — One-dimensional H+K band spectra for the two components of J0900+2234.The upper panel is the spectrum of component A 
and the lower panel is the spectra of component B. The dashed line is Icr error of the spectrum, and the dash dot line represent the zero 
flux. 
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Fig. 4. — g, r, i, J, H and Ks-haxid AB magnitudes of the knot A (the upper panel) and the knot B (the bottom panel) and their best-fit 
model spectra.The best fit age is 180 Myrs for both knots and the lensed corrected stellar mass is 1.9 x 10^°Mq. 
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TABLE 1 

LUCIFER Observation Log of J0900+2234. 



Datc(UT) 


mode 


Camera 


Filter 


Slit PA (degree) 


Exp. Timc(s) 


Jan. 06 2010 


imaging 


N3.15 


J 




10 X 60 








H 




10 X 60 








Ks 




10 X 60 


Feb. 13 2010 


spectroscopy 


N1.8 


H+K 


86.57 


24 X 300 



TABLE 2 

The best-fit profile parameters for the lens galaxies and lensed gomponents. 



Component 


Model 


Effective Radius (-Re(")) 


Axis ratios ib/a) 


Position angle (9(°)) 


01 


de Vaucoulcurs 


4.19 ±0.08 


0.70 ±0.01 


-60.63 ± 1.00 


G2 


de Vaucouleurs 


1.52 ±0.03 


0.73 ±0.01 


7.16 ± 1.47 


A 


exponential disk 


0.45 ± 0.04 


0.36 ±0.03 


-17.84 ±3.22 


B 


exponential disk 


0.55 ± 0.04 


0.58 ± 0.04 


-15.70 ±4.93 


V 


exponential disk 


0.45 ± 0.02 


0.71 ± 0.03 


48.90 ± 5.48 


D 


exponential disk 


2.39 ±0.27 


0.11 ±0.01 


46.52 ±0.95 


Al 


exponential disk 


0.63 ± 0.05 


0.34 ± 0.02 


12.40 ± 2.38 



TABLE 3 

Photometry of the LUCIFER images from GALFIT. 



Component 


Jab 


Hab 


KSAB 


Gl 


17.34 ±0.01 


16.98 ±0.01 


16.68 ± 0.03 


G2 


18.43 ± 0.01 


18.06 ± 0.01 


17.69 ± 0.02 


A 


20.99 ± 0.03 


20.87 ± 0.04 


20.87 ±0.06 


B 


20.63 ± 0.02 


20.33 ± 0.04 


20.30 ± 0.04 


C 


20.45 ± 0.03 


20.44 ± 0.04 


19.97 ±0.03 


D 


20.79 ±0.07 


21.08 ±0.18 


20.23 ± 0.09 


Al 


21.81 ±0.09 


21.26 ±0.11 


20.47 ± 0.04 



TABLE 4 

Results of emission line measurements. 



Components 


Lines 


wavelength 


flux 


FWHM 






A 


10"^'' ergs s~^ cm~2 


A 


A 


A4863 


14738.9 


4.85 ± 1.78 


18.2 ±8.2 


B 


H/3 A4863 


14738.1 


15.08 ± 1.36 


18.7 ± 1.4 


A 


[O III] A4959 


15043.0 


14.15 ± 2.99 


22.5 ±8.1 


B 


[O III] A4959 


15037.4 


29.98 ±3.93 


23.1 ±2.6 


A 


[O III] A5007 


15184.6 


37.24 ± 1.78 


22.4 ± 1.1 


B 


[O III] A5007 


15182.7 


93.81 ±4.07 


26.4 ± 1.1 


A 


Hq A656:? 


19902.8 


37.01 ±0.79 


19.0 ±0.4 


B 


Hq A6563 


19899.5 


84.68 ±3.21 


21.8 ±0.6 


A 


[N II] A6584 


19965.1 


1.56 ± 0.35 


14.6 ±4.3 


B 


[N II] A6584 


19961.9 


5.66 ± 0.88 


22.3 ±3.8 


A 


[S II] A6717 


20367.5 


3.56 ± 0.50 


22.0 ±3.9 


B 


[S 11] A6717 


20365.1 


5.11 ± 1.01 


19.9 ±5.6 


A 


[S II] A6732 


20405.6 


4.06 ± 0.98 


17.3 ± 5.3 


B 


[S II] A6732 


20410.4 


6.06 ± 1.87 


20.6 ±8.6 



TABLE 5 

Summary of the physical properties of J0900+2234. 



Knots 


age 


Es(B-V) 




Zk2 


•Z03N2 




SFRhc" 


SFRuv" 


loglM^tellar)" 


log(Mgasr 






Myr 






Zq 


Zq 


cm^"^ 


Mq yr-1 


Me yr-i 


Mq 


Me 


Me 


A 


180 


0.07 


0.84 ±0.31 


0.27 + 0.13 


0.21 ± 0.08 


I029ti];];f 


365 ± 69 


203 ± 38 


10.28 


10.71 


10.76 


13 


180 


0.20 


0.51) ± 0.08 


O.-'i") = 0.15 


0.2() = 0.09 


il(i(j_;j',|j" 













^This value is derived from average of A and B components. 



